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Free energy of an SU(2) monopole-antimonopole pair * 
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We induce an external Z2 monopole-antimonopole pair in an SU(2) lattice gauge system and measure its free 
energy as a way to probe the vacuum structure. We discuss the motivation and computational methodology of 
the investigation and illustrate our preliminary results. 
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It is well known that some Higgs theories with 
non-Abelian gauge group admit stable monopolc 
solutions H,||]. With large gauge groups, as m 
grand unified theories, the residual unbroken 
gauge group can be non-Abelian. It is then in- 
teresting to determine the properties of the inter- 
action induced among monopoles, or, as we will 
consider in this paper, beteween a monopole and 
an antimonopole, by the quantum fluctuations of 
the unbroken group. Beyond the relevance that 
such interaction may have for the original theory, 
it can shed light on the low energy properties of 
the residual gauge theory itself. From the point 
of view of the unbroken theory, to a very good 
approximation the monopoles act as static point 
sources. The way to incorporate such sources in 
an SU(N) lattice gauge theory was spelled out 
in Refs. ||(|, which built on earlier results es- 
tablished in a seminal paper by 't Hooft [0 and 
in Refs. We will follow closely the treat- 

ment of Ref. j|. In a three dimensional theory, 
a monopole-antimonopole pair can be introduced 
within two cubes of the lattice by the replacement 
PTtUp -> P'TtUp = z„/?Tr[/ P for all the terms 
in the action corresponding to plaquettes trans- 
versed by a path joining the centers of the two 
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cubes. z n stands here for a non-trivial element 
of the center of the group SU (N) i.e. for one of 
the N th roots of unity: z n = ex-p(2mn / N) with 
n = 1 ... N — 1. In the case of SU(2) which we 
will study here the only non-trivial element of the 
center of the group is -1 and monopole and anti- 
monopole coincide. By the redefinition U — > z n U 
of an appropriate set of the link variables, the 
path joining monopole and antimonopole can be 
deformed at will. The path is therefore unphysi- 
cal and carries no free energy per se. But its end 
points cannot be moved without an accompany- 
ing change of free energy. A static monopole- 
antimonopole pair is introduced in the four di- 
mensional theory by simply replicating the above 
construction in all time slices. Now the plaquettes 
with modified coupling are those transversed by 
a sheet joining the world lines of the monopoles. 
Again the sheet can be deformed at will. The in- 
sertion of a monopole-antimonopole pair can be 
rcintrepeted in terms of the electric flux operator 
introduced in Ref. 0] and several investigations 
have been devoted to the study of such operator 
in various contexts (see for example Refs. PJlOfl), 
but, to the best of our knowledge, no direct cal- 
culation of the free energy of a Z2 monopole- 
antimonopole pair has ever been attempted. If 
the vacuum is characterized by a monopole con- 
densate, one would expect the behaviour of the 
free energy as function of separation to exhibit 
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screening, whereas a 1/r behavior or a linearly 
rising behaviour would characterize a Coulomb 
phase or a phase with condensation of electric 
charges, respectively. It is also noteworthy that 
the monopole and the antimonopole form two an- 
chors for a center vortex. Recent investigations 
(cfr. Refs. have emphasized the role that 

such vortices play in confinement. The calcula- 
tion which we present here can be reinterpreted as 
the calculation of the cost in free energy to create 
a center vortex spanning a certain distance within 
the lattice. If such excess free energy quickly satu- 
rates (screening) , then the vacuum should indeed 
exhibit a condensate of center vortices. 

The numerical calculation of a free energy is 
notoriously difficult. We have been able to ob- 
tain reasonably accurate results with acceptable 
amounts of CP time by combining a Monte Carlo 
simulation with the multihistogram method |l5| . 
We consider a modified SU (2) lattice gauge the- 
ory with Wilson action, defined over a N x x 
N y x N z x N t hypercubical lattice with peri- 
odic boundary conditions. The modification con- 
sists in the fact that, for all the x — y plaque- 
ttes P' having a lower vertex with coordinates 
x = 0, y = 0, < z < d, < t < Nt, the coupling 
constant (3 is replaced with f3' . These are the pla- 
quettes that cross the sheet joining the worldlines 
of the monopole and antimonopole at separation 
r = da (a being the lattice spacing). We denote 
the partition function of this system by Z(f3',f3). 
We are interested in the free energy 
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If we perform a simulation with (3' set to a certain 
value (3i and record in a histogram the frequency 
rii(E) of occurrences of a certain value of E, we 
will find 

^From p{E) (as yet unknown) we get 
Z{Pi,P)= [ dEp(E)e^ E/2 (4) 



The procedure we have followed consists in per- 
forming a number Nh of simulations covering 
the interval —(3 < /3' < /3 with a sufficiently 
small step to guarantee a good overlap of the his- 
tograms. We take then 
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1 and iterate: the values of 
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z[ k ^ at convergence are proportional to the cor- 
responding Z(j3i, (3). As a technical improvement 
in the implementation of the histogram method, 
applicable to the case where the measured vari- 
able covers a continous range, we have allocated 
the values of E to the four neigboring end-points 
of the histogram intervals with the weights of a 
cubic interpolation. This procedure reduces sub- 
stantially the total number of histogram subdivi- 
sions one must use to obtain accurate results. 

We illustrate here the results we have obtained 
with (3 = 2.6, N x = N y = 20, N z = 40 and the 
two time extents Nt = 16 and Nt = 6, which place 
the system in the confined and deconfined phases 
respectively. We have used a combined multi-hit 
Metropolis overrelaxation algorithm, with 5000 
equilibrating iterations and 4000 to 20000 mea- 
surements separated by 50 iterations. The mea- 
surements themselves have been performed by av- 
eraging over 384 upgrading steps of the links in 
the plaquettes P' as a variance reduction tech- 
nique. In Figure 1 we show the histograms for a 
definite separation of the monopole- antimonopole 
pair. In Figure 2 we show our results for the free 
energy of the pair. The calculation is computer 
intensive, because one must perform separate cal- 
culations for all separations of the pair and for all 
the intermediate values of (3' . However it is quite 
feasible with present day computer resources. We 
have written a Fortan 90 code which, paying some 
attention to the distribution of the data but with- 
out resorting to any special programming tricks 
like coding critical subroutines in assembler, runs 
at approx. 40% of peak speed on the SGI-Cray 
Origin 2000, with very satisfactory scaling. With 
this performance, the cost of the data presented 
in this paper is of the order of a few thousand pro- 
cessor hours, which is rather modest by today's 
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Figure 1. The overlap of histograms at (3 = 2.6 
and different /?'. The monopole separation is 2a 
and N t = 6. 
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monopole pair separation 

Figure 2. Free energy of a monopole- 
antimonopole pair at different separations for f3 — 
2.6 and N x x N y x N z = 20 x 20 x 40. The cir- 
cles are with Nt — 16 in the confined phase, the 
diamonds with Nt = 6 in the deconfined phase. 



standards of supercomputing. 

The results in Figure 2 show that the inter- 
action of the pair is screened both in the con- 
fined and deconfined phases. The lines in the 
figure correspond to exponential fits exp(— d/l) 
with I = 0.8300 and I = 0.7828 for N t = 16 and 
Nt = 6, respectively. We have also computed 
the free energy of a single monopole adopting free 
boundary conditions for the z = 0, N z boundaries 
of the lattice (we maintained periodic boundary 
conditions in all other directions) and the results 
are in agreement with the free energy of the pair 
for large separation. Our investigation is still 
in progress. We plan to repeat the calculation 
with a smaller value of f3 to verify scaling and to 
study the behavior of the free energy of a single 
monopole as one goes across the deconfining tran- 
sition. It would also be interesting to extend the 
calculation to other systems, especially to models 
which are expected to possess a varied structure 
of electric and magnetic confinement phases. 
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